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Abstract. In this paper we focus on the pathwise stability of mild solu- 
tions for a class of stochastic partial differential equations which are driven by 
switching-diffusion processes with jumps. In comparison to the existing liter- 
ature, we show that: (i) the criterion to guarantee pathwise stability does not 
rely on the moment stability of the system; (ii) the sample Lyapunov exponent 
obtained is generally smaller than that of the counterpart driven by a Wiener 
process; (iii) due to the Markovian switching the overall system can become 
pathwise exponentially stable although some subsystems are not stable. 



1. Introduction 

Stochastic partial differential equations (SPDEs) have been widely used to model 
phenomena arising in many branches of science such as ecology, economics, mechan- 
ics, biology and chemistry, e.g., Applebaum [3J, Chow [9], Da Prato and Zabczyk 
[TO] . Liu [16], Peszat and Zabczyk [21], and Woyczyhski [25] , Recently, hybrid 
systems, in which continuous dynamics are intertwined with discrete events, have 
also been used to model many such systems. One of the distinct features of hybrid 
systems is that the underlying dynamics are subject to changes with respect to 
certain configurations. For example, consider a one-dimensional rod of length it 
whose ends are maintained at 0° and whose sides are insulated. Assume that there 
is an exothermic reaction taking place inside the rod with heat being produced 
proportionally to the temperature. The temperature u in the rod may be modelled 
by (see, e.g., [11]) 



(1.1) 



du - d u -^CU, t > 0, X € (0,7T), 



dt dx 

u(t,0) = u(t,ic) = 0, 
u(0,x) = u (x), 



where u = u(t, x) and c is a constant dependent on the rate of reaction. The 
system (jl.ip will switch from one mode to another in a random way when it expe- 
riences abrupt changes in its structure and parameters caused by phenomena such 
as component failures or repairs, changing subsystem interconnections, or abrupt 
environmental disturbances. A hybrid system driven by a continuous-time Markov 
chain can be applied to describe such a situation. The system ([l.ip under regime 
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switching could be described by the following stochastic model 

'§ = ^ + c(r(t))u, t>0, xG(0,n), 
< u(t,0) = u(t,ir) = 0, 
u(0, x) = u (x), r(0) = r , 

where r(t) is a right-continuous Markov chain with finite state space § and c : 
§ — > R. As a second example, Li et al. [H] discussed stochastic age-dependent 
population equation with Markovian switching, 

-/i(r(t), a)P + f(r(t), P) + g( r (t),P)^, (t, a) e Q, 
P (a), r(0)=r , o6[0,4 
f*P(t,a)P(t,a)da, t e [0,T], 

where T > 0,A > 0,Q := (0,T) x (0,A),P = P{t 7 a) is the population density 
of age a at time t, r(t) is a right-continuous Markov chain, fi(r(t),a) denotes the 
mortality rate of age a at time t, and f(r(t),P) denotes the effects of the external 
environment on the population system. For more quantitative analysis of SPDEs 
with Markovian regime switching, we refer to Anabtawi and Ladde ]T|, Anabtawi 
and Sathananthan [2 ], Luo and Liu |17j . and the references therein. For the finite- 
dimensional case, we refer to the monographs of Mao and Yuan [18] and Yin and 
ZhugSJ. 

Non-Gaussian random processes also play an important role in modelling sto- 
chastic dynamical systems, e.g., Applebaum [3], 0ksendal and Sulem [21], and 
Peszat and Zabczyk [24]. Typical examples of non-Gaussian stochastic processes 
are Levy processes and processes arising from Poisson random measures. The 
monograph [23] describes a number of phenomena from fluid mechanics, solid state 
physics, polymer chemistry, economic science, etc., which can be modelled using 
non- Gaussian Levy processes. 

Moreover, one of the most important and interesting problems in the analysis of 
SPDEs is their stability. Stability issues for SPDEs with Wiener noise are by now 
classical, see, e.g., [13l HH [17], but comparable theories driven by jump noise 
are not yet fully developed. Recent years have witnessed a growing interest in this 
area: In [6] we investigated the asymptotic stability in distribution of mild solutions 
for delay equations using Lyapunov functions and the Yorsida approximation; by 
energy inequality approach, moment stability and sample path stability for varia- 
tional solutions were discussed in [12] . and almost sure exponential stability was 
studied in |17| provided that the mild solution is moment exponentially stable. For 
the finite-dimensional case, Applebaum and Siakalli [4 provide sufficient conditions 
under which the solutions to stochastic differential equations (SDEs) driven by Levy 
noise are stable in probability, almost surely and moment exponentially stable. As 
a sequel, Applebaum and Siakalli [5] made some first steps in the stochastic stabi- 
lization problems where the noise source is a Levy noise, i.e., a Brownian motion 
and an independent Poisson random measure. 

Most of the previous literature does not discuss the sample Lyapunov exponents 
which are given explicitly by the parameters arising from the jump-diffusion coef- 
ficients or stationary probability distribution of Markovian chain. In this paper, 
motivated by the previous references, we shall study the pathwise stability of mild 



dP I OP _ 
dt " r da 

P(0,a) = 
P(t,0) = 
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solutions for a class of SPDEs of the form 

dX{t) = [AX(t) + F(t, X(t), r{t))]dt + G(t, X(t),r(t))dW(t) 

(1.2) f 

+ <Z>(t,X(t-),r{t),u)N(dt,du) 
Jz 

on t > with the initial data X(Q) = xq £ H and r(0) = tq e S, where X(t~) := 
lim s -ft X(s). More detailed information on the parameters in Eq. (|1.2[) will be given 
in Section 2. 

In comparison to the existing literature for the almost sure exponential stability 
of solutions to SPDEs, our main result (Theorem 3.1), has the following advantages: 
(i) the criterion established does not rely on the moment exponential stability of 
the system; (ii) the sample Lyapunov exponent is generally smaller than that of 
the counterpart driven by a Wiener process; (iii) due to the Markovian switch- 
ing the overall system could become pathwise exponentially stable, although some 
subsystems are not stable. 

Our approach is based on the Yosida approximation and a classical Lyapunov 
function argument. Since mild solutions do not necessarily have stochastic differen- 
tials, one cannot apply the Ito formula directly. To overcome this problem, we first 
apply the Ito formula to an approximating equation, and then investigate the sta- 
bility properties of the mild solutions. This approach is dependent on a maximal 
inequality (Burkholdcr-Davis-Gundy inequality) for stochastic convolutions with 
jumps. 



2. Preliminaries 

Let (H, (•, -)h, II • \\h) be a Hilbert space and W(t) a cylindrical Wiener pro- 
cess on H defined on some filtered probability space (fijJ 7 , P) equipped with a 
filtration {J r t}t>o satisfying the usual conditions. Let Z be a vector space en- 
dowed with a norm | • |, £>(Z) is the Borel tr-algebra on Z, and \{dx) a a- finite 
measure defined on £>(Z). Let p — (p(t)),t £ D p , be a stationary J^-Poisson 
point process on Z with characteristic measure A. Denote by N(dt, du) the Pois- 
son counting measure associated with p, i.e., N(t,Y) = J2 s eD p s<t ^(p( s )) f° r 
Y £ B(Z). Let N(dt,du) := N(dt,du) — dtX(du) be the compensated Poisson 
measure. Let m be some positive integer and {r(t),t £ K+} a right continuous ir- 
reducible Markov chain on the probability space J 7 , P} taking values in a finite 
state space S := {1, 2, m}, with generator T = (7y)mxm given by 



P(r(t + A) =j\r(t) =i) = 



7yA + o(A), if i^j, 

l+ 744 A + o(A), if i = j, 



where A > and 7^ > is the transition rate from i to j, if i ^ j; while 
7« = ~ Sj^i lij- F° r a mapping £ : § — > R, we write Q := We further assume 
that the Wiener process W(t), Poisson process N(t, •) and Markov chain r{t) are 
independent. For T > and p > 0, let -D([0, T];H) be the family of all cadlag paths 
from [0,T] into H with the uniform norm and L' } := L p (Vt, J 7 , P; D([0, T]; H)) := 
{X:Esu Po < t < T ||X(t)||^<oo}. 

Remark 2.1. Since we have assume that the Markov chain r(t) is irreducible, it has 
a unique stationary probability distribution ir := (tti, 7T2, . . . ,Trm) £ M lxm which 
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can be determined by solving 

I a 

ttT = subject to iij = 1 and ttj > 0, Vj G S. 

3=1 

Throughout the paper we will assume: 

(HI) A : T>(A) — > H is the infinitesimal generator of a pseudo-contraction Co- 
semigroup S(t) on H, that is, 

(2.1) \\S(t)\\ < e Qi , for some a > 0. 

(H2) F : [0,T]xiJx§^ H, G : [0, T] x H x § ->• £ 2 (-H, the family of Hilbert- 
Schmidt operators from into itself, and $ : [0, T] x _ff x Z x § — > H , satisfy 
Lipschitz and linear growth conditions, i.e., there exist positive constants 
L, L such that for arbitrary x,y £ H,i 6 § 

\\F{t,x,i)-F(t,y,i)\\ 2 H + \\G{t,x,i)-G(t,y,i)\\ 2 HS < L\\x -yf H , 
\\F(t, x, i)f H + \\G(t,x,i)\\ 2 HS< L(l + \\x\\ 2 H ), 

/ \\$(t,x,i,u) - $(t,y,i,u)\\ 2 H \(du) <L\\x-y\\ 2 H , 

f \mt,x,i,u)\\ 2 H \(du) <L(l+\\x\\ 2 H ). 

Moreover, we will need some slightly stronger conditions on $: 
(H3) For p > 2, there exists L > such that 

(2.2) [ \\<f>{t,x,i,u)-$(t,y,i,u)\\ p H \(du) <L\\x-yf H , x,yeH,ie§, 
Ji 

and <f>(t,Q,i,u) = 0. 

(H3 1 ) There exists an L > such that for arbitrary x,y £ H,i G S,u G Z and 
p > 2 

||$(i, a;, i, u) — $>(t,y, i,u)\\ H < L\\x — and / \u\ p X(du) < oo. 

Clearly (H3') implies (#3) and ||$(£, a;, i, < L||a;||jj|u| whenever $(£, 0, i, u) — 
0. 

In this paper we are mainly concerned with mild solutions to (jl.2p . For the 
notion of mild solutions, we can refer to, e.g., 19, Definition 2.1], for SPDEs driven 
by multiplicative Poisson noise, and, with an obvious extension, we can define such 
solutions for the case of switching-diffusion processes with jumps. 

Lemma 2.2. Under conditions (HI) and (H2), Eq. (|1.2j) admits a unique mild 
solution X(t,x ,r ),t G [0,T], in L 2 (Q, J 7 , P; H). 

Proof. Since the proof of Lemma 12.21 can be done by combining classical Banach 
fixed-point theorem argument [24J Theorem 9.29, pl64] with stopping time tech- 
niques [18[ Theorem 3.13, p89], we will only give a sketch of the argument. Recall 
that almost every sample path of r(-) is a right-continuous step function with a 
finite number of sample jumps on [0, T\. So for almost every uj G there is a finite 
k = k(uj) such that = To < T\ < ■ • ■ < Tu > T and 

r(t) = r(r k ) on r fe < t < T k+1 for Vfc > 0. 
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We first consider Eq. (|1.2p on t G [0, ri] which becomes 

dX(i) = [AX"(i) + F(<, X(t), r )]dt + G(t, X(t),r )dW(t) 
$(t, X(i~), r 0) u)N(dt, du) 



(2.3) 



with the initial data X(0) = zo € H and r(0) = r e S. By [211 Theorem 9.29, 
pl64] Eq. (|2.3[) admits a unique mild solution X(t,xo,ro) for t 6 [0, n] which 
belongs to L 2 (fl, J 7 , P; H). Next consider Eq. (|3.ip on i 6 [n,^] which becomes 

dX(t) = [AX(t) + F(t, X(t), r(Ti))]tft + G(t, X(*), r(n))(flf(t) 

(2-4) 

' $(t, r(T 1 ),u).N'( £ Mu) 

with the initial data X(ti, x n ,r ) and r(ri). Again by [Ml Theorem 9.29, pl64] 
Eq. (|2.4p admits a unique mild solution X(t, X(tx, Xo, T"o), t(ti)) for f G [ti,t 2 ] 
which belongs to L 2 (Vt, J", P; H). The proof can then be completed by repeating 
this procedure. □ 



Let U : R+ x H x § -> R+ be a C^-function. For t > 0, x E V(A) and i 6 S, 
define an operator 

CU(t,x,i) := U t (t,x,i) + (Ax + F(t,x,i),U x (t,x,i)) H 

m ^ 

+ ^2 ^v U ( t ' x ' i) + ^^ce(U xx (t, x, i)G(t, x, i)G* (t, x, i)) 
i=i 

+ / x + $(t, x, i, u), i) - U(t, x, i) - (U x (t, x, i), $(t, x, i, u)) H ]\(du). 
Jz 

Similarly to [H Theorem 1.45, p48] and [H Theorem D.2, p392], for a strong 
solution X(t) of Eq. (jl.2[) . we have the following ltd formula: 

(2.5) 

U(t,X(t),r(t)) 

= (7(0, xq, r ) + / £C/(s, X(s), r(s))ds 
Jo 

+ f (U x (s,X(s),r(s)),G(s,X(s),r(s))dW(s)) H 
Jo 

+ f f [U(s, X(s-) + X(s-), r(s), u),r(s)) - U(s, r( S ))]iV(ds, du) 

Jo Jz 

+ / /[C/( S ,X( S -),r + ^(r( S ),f))-C/( S ,X( S -),r( S )) A1 (d S ,^), 

where fi(ds,d£) is a Poisson random measure with intensity ds x m(d£), in which 
m is the Lebesgue measure on R. For more details on the function h and the 
martingale measure /j,(ds,d£), see, e.g., JT5J p46-48]. 

Since mild solutions do not necessarily have stochastic differentials, one cannot 
apply the Ito formula directly. Instead, we introduce the Yosida-approximation 
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system 

dX l (t) =[AX l (t) + R(l)F(t,X l (t),r(t))}dt + R(l)G(t,X l (t),r(t))dW(t) 

(2.6) 



R(l)<f>(t, X l (t~), r(t), u)N(dt, du) 



'z 

with the initial data X l (0) = R(1)xq and r(0) = r . Here I e p(A), the resolvent 
set of A, and R(l) := IR(1,A) := l(lld — A) -1 , where Id is the identity operator 
from H into itself. 

Theorem 2.3. Under conditions (HI) — (H3), Eq. (|2.6[) has a unique strong 
solution X l (t) in > 2. Moreover, X ( (i) converges to the mild solution X(t) of 
Eq. (|L2|) in L p as I ^ oo, i.e., 

EmE sup ) = 0. 



. te[o,r] 



In particular, there exists a subsequence, still denoted by X l (t), such that X l (t) — > 
-X"(t) almost surely as I — > oo uniformly in [0,T]. 

Proof. This can be proven by following the argument in |17[ Proposition 2.4] and 
using the result of Lemma 12.41 below. □ 







sup 


III 


0<t<T 





Lemma 2.4. ( [7J Theorem 4.4] or [TH Proposition 3.3]) Assume that $ : fl x 
K + x Z — > H, is a progressively measurable process, and for p > 2 

(2.7) E f [ \\<$>(s,u)\\ p H \(du)ds < oo. 

Jo Jz 

If S(t) 7 t e [0,T], is a pseudo-contraction Co-semigroup such that (|2.1j) . then 
(2.8) 

r* /■ p 1 /" T /" 

S*(t-s)$(s,u)A>(ds,du) <C P E \\$(s,u)\\ p H \(du)ds, 
: h Jo J j, 

where C p is a positive constant dependent on p, a, T. 

3. A Criterion for Sample Path Stability 

In this section we give a criterion for sample path stability of mild solutions to 
Eq. (jL2l . 

Theorem 3.1. Let (HI), (H2) and (Hi 1 ) hold. Assume that the solution of Eq. 
(|1.2j) is such that X(t) 7^ a.s. for all t > and i G 8 provided xo ^ 
a.s. for [7 £ C 1,2 (K+ x H x §;R + ), assume further that there exist constants 
c 2 > ci > 0,c 3 > 0,p > 0, ai, pi S R, ft, 5i > smc/i that for (t, x, i) e R + xV(A)xS> 



(i) ci||i||* < U(t,x,i) < c 2 \\x\\ p H , \\U x (t,x,i)\\ H \\x\\ H + \\U xx (t,x,i)\\\\x\\ 2 H < 

csh\\ p H ; 

(ii) £U(t,x,i) < ctiU(t,x,i); 

(hi) SU(t,x,i) := ||G*(i,a;,i)[4(i,a:,i)|& > P>U 2 (t,x,i); 
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(iv) For #(t, x, i, j) := U(t, x, j)/U(t, x, i), 

m 

j'=i 

(v) For A(t, x, i, u) := C7 (t, a; + x, i, u), i)/U(t, x, i), 
[lnA(t, x, i, u) — A(t, x, i, u) + 1] X(du) := J(t, x, i) < — 5j, 

and for some e £ (0, s] 

£ := limsup — H(t) < oo, 

t— >oo t 

where 

E(i) := / / [(In A(s,x,i,w)) 2 + A e (s, x, i, ■u)]A(du)ds; 
Jo ./Z 

(vi) For T(t,x,rQ,i,£) :— U(t,x,ro + h(i,£))/U(t,x,i) and some e £ (0, i] 

?7 := limsup — < oo, 

t— >00 i 

w/iere 

ct 



n(t) := / / [(In T(s, x, r , i, ^)) 2 + T e (s, x, r , i, £)]m(d£)d 



s. 



/o Jz 

Tften i/ie miH solution of Eq. (11.21) /ias the property 

1 1 m (\ 

limsup -ln(||A(t)|| H ) < — V tt 2 -/3, + (5, - on - 
t-Voc t P~[ \ 2 

In particular, the mild solution of Eq. (|1.2p is almost surely exponentially stable 
provided that 

i=l ^ ' 

Remark 3.2. By the fundamental inequality 

ln(l + x) < x for x > 0, 

it is easy to observe that the first assumption in (v) is reasonable. On the other 
hand, for x £ H and i, j £ S, (vi) is also true provided that there exist constants 
C4, C5 > such that 

c 4 < U(t,x,j)/U(t,x,i) < c s . 

There are many functions possessing this property, e.g., for x E H,i £ S, U(t, x, i) = 
<7i||x||jy with Ui > 0. For this Lyapunov function, we have ^(t, x,i, j) = (Tj/cri and 
then (iv) also holds. 

Proof of Theorem 13. li As we mentioned previously, one cannot apply the Ito 
formula to the mild solutions directly, we first apply the Ito formula to the ap- 
proximation equation (|2.6|) . then use Theorem [273] to investigate the stability of the 
mild solutions. More precisely, applying the Ito formula (|2.5|) to In U(t, x, i) with 
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respect to X l (t),t > 0, where X l (t) denotes the strong solution of Eq. (|2.6j) . we 
obtain 



In U(t, X l (t), r(f)) = In £7(0, R(l)x , r ) 

rt m 



where 



1 CU(s,X l (s),r(s)) 



ds 



o U(s,Xi(s),r(s)) 

/ I>r«:i [ln*(s,X l (s),i,j)-*(s,X l (s),i,j)]ds 
Jo j=1 

i ' In ( U ^ X ^ tg^gfe *' ( '~ >' r(5) ' r(5)) ) N(ds, du) 



U(s,X l (s-),r(s)) 



\nr(t,X l ( S -),r ,r{s),£)[i(ds,d£) 



+ J 1 (t,l) + J 2 (t,l) + J 3 (t,l), 



,h{t,l) 



MM) 



trace 



U xx (s, X l (s),r(s)) U x {s, X l (s), r(s)) ® U x (s, X l (s), r(s)) 



U(s,Xi(s),r(s)) 



U 2 (s,Xi(s),r(s)) 



R(l)G(s, X l (s), r(s))(R(l)G(s, X l (s),r(s))Y 



ds 



t trace 



U xx (s, X l (s), r(s))G(s, X l (s), r(s))(G(s, X 1 (s), r(s))Y 



U(s,Xl(s),r(s)) 
(U x (s,X l (s),r(s)),(R(l) - I)F(s,X l (s),r(s))) H 



ds 



U(s,Xl(s),r(s)) 



ds. 



In 



' U(s, X l (s) + RjlMs, X l (s), r(s), u),r(s)) 

U(s,Xl(s),r(s)) 
(U x ( S ,X l (s),r(s)),(R(l) - lMs,X l (s),r( S ),u),r(s)) 

U(s,Xl(s),r(s)) 
U(s,X l (s) + *(s,X l (8),r(s),u),r{8)) 
U(s,X\s),r(s)) 



X(du)ds. 



By Theorem 12.31 (H3'), (i) and the Dominated Convergence Theorem, we have 
almost surely 



1 [' eU(s,X( 3 ),r(s)) j 

km Ji (*,/) = -- / — — ds, hm. J 2 {t,l)=0, 

z-s-oo I J U (s, X{s),r{s)) i-s-oc 

lim J 3 (t,I)= / J(a,X(8),r(s))ds. 
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If we let / — > oo, then, 



\nU{t,X(t),r(t)) =\nU{0,x o ,r a ) + / J{s, X(s),r(s))d 



+ 



CU(s,X(s),r(s)) _1 f* eU(s,X(s),r(s)) 
„ U(s,X(s),r(s)) 3 2j U*(s,X(s),r(s)) 
4 (U x (s, X(s), r(s)),g(s, X(s), r(s))dW(s)) H 



ds 



U(s,X(s),r(s)) 
/ X)7r(.)j [ln*(s,X(s),i,j) - *(«,X(«),i,i)]dfl 

^0 •_! 



j = l 



/o Jz 



In A(s, X(s~), r(s), u)N(ds, du) 
In T(s, X(s"), r , r(s), %(ds, 



By the exponential martingale inequality with jumps [3j Theorem 5.2.9, p291], for 
any T,9,v>0 



uj : sup / 

0<t<T L Jo 



(^( S ,X( S ),r( S )),. 9 ( S ,X( S ),r( S ))^( S )) g ^ /* eC/(,s,X(,s),r( S )) 
[/"(«, X(*),r(a)) 



2 7 C7 2 ( S ,X( S ),r(s)) 



+ / / lnA(s,X(s _ ),r(s),it)iV'(ds,(iu) 







< e 



Jz 

-8u 



A e (s,X(s),r{s),u) - 1 - 9\nA{s,X(s),r{s),u) X(du)ds > i/j 



and, if we denote T(t,£) := Y(t, ie, ro,i,£), then 



w : sup 

L 0<t<T 



lnf (s,£)ii(ds,de) 



'(s,£)-l-01nf(s,^) m{d£)ds > ^| 



< e 



-01/ 



Taking T = n, v= 29~ x Inn, n = 1,2,..., for £ (0, ^] and applying the Borel- 
Cantelli Lemma, we see that there exists an f?o C f2 with P(ilo) = 1 such that for 
any w G Jlo there exists an integer hq — no(uj) > if n > no 

1 {U x (s, X(s),r(s)),g(s, X(s),r(s))dW(s)) H " 



U(s,X(s),r(s)) 



In A(s, X(s~), r(s),u)N(ds, du) 



o Ji 



< 2<T 1 lnn+ — 



9 2 r eU(s,X(s),r(s)) 



2 J C/ 2 (s,X( s ),r( S )) 



ds 



A e {s,X(s),r(s),u) - 1 -6\nA(s,X(s),r(s),u) X(du)ds 



and 
ft 

/o 



lnf (s,l)ii{ds,dl) < 2<r 1 hr, 



f e (s,£) - 1 -6» lnf (s,£) m{dt)ds 
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for < t < n. Hence for cjG£lo,0<i<n and n > no 

lnU(t,X(t),r(t)) < lnE/(0,x o ,r o ) + 4<r 1 lnn + I J{s,X(s),r(s))ds 

Jo 

'* QU(s,X(s),r(s)) 



* CU(s,X(s),r(s)) 1 
as 



o U(s,X(s),r(s)) 



2 J U 2 (s,X(s),r(s)) 



ds 



„t m 

[ln^(s,X(s),i,j)-^(s,X(s),i,j)]ds 



1 

+ e 
i 



Ji 
t 



JR 



A e (s,X(s),r(s),u) - 1 - 0ln A(s,X(s),r(s),u)\\(du)ds 
f (s,£) - 1 - eint (s,l)]m(dl)ds. 



This, together with (ii)-(v), yields that for uj e SIq , < t < n and n > no 
hiU(t,X(t),r(t)) < ln[/(0,x o ,r o ) +4(T 1 hin 

9 



+ 



f [a(r(s))-- 
Jo 



-P(r(s))-6(r(s))+p(r(s))]ds 



(3.1) 



T(s,X(s),y,r(s),e)X(dy)ds 



io Ji 

rt 



where 
and 



e(s,£,6)m{dt)ds 
:= h(t) + I 2 {t) + I 3 (t) + I 4 (t), 

T(s,x,i,u,0) :— A e (s,x,i,u) — 1 — 9 In A(s, x, i, u) 



@{s,i,6) ■=T e {s,t)-l-9\nT{s^). 
Using Taylor's series expansion, for 9 e (0, ^] we see that 

9 2 

A 8 (s,x,i,w,) = 1 + #lnA(s,x,i,u) + — (In A(t, x, i, u)) 2 A^ (t, x, i, u) , 
where £ lies between and 9. Hence 

h{t)= 6 - f I (lnA(s,X(s),r( S ),u)) 2 AH S ,X( S ),r(s),u)\(du)ds 
z Jo Jz 

= 9 - [ f (lnA(s,X(s),r(s),u)) 2 At(s,X(s),r(s),u)\(du)ds 
2 Jo J0<A<1 

+ e - [ f (lnA(s,X(s),r(s),u)) 2 A i (s,X(s),r(s),u)X(du)ds. 
2 Jo Ja>i 

Noting that, for < £ < § , A« < 1 if < A < 1, A« < A 5 if A > 1, and recalling 
the inequality 

\nx < -(x* — 1) for x > 1, 



we obtain 



77 

Jo Jz 



(In A(s, A(s), r(s), u)) 2 + ^A e (s, A(s), r(s), u)] X(du)ds. 
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Similarly, 



h(t) < - 



16 i 

(lnT(s^)) 2 + —T e (s,£) m{d£)ds. 



Hence, by (i) for u 6 fio> n — I < t < n and n > uq + 1 

- H\\X{t)\\ H ) < + - [in £7(0, xo, r ) + 4^ Inn + r -^ I (H(i) + U(t)) 

t pt ptl (e A e) J 



+ / Wr(s)) - 



I- a 



p(r(s)) ~ 6(r(s)) + p{r(s))]ds 



Taking into account the ergodic property of Markovian chains, e.g. 
3.8.1, pl26], and combining (v) with (vi), we have almost surely 



Theorem 



nmsup~m(||X(t)|| H ) < - 

t^oc t p 



(eAi) : 



:(C + »») + $>< 



1 



-f3. t -Si + p, 



and the conclusion follows by the arbitrariness of 9. 
Now we give an example to demonstrate Theorem 13.1 



Example 3.3. Let r(t) be a right-continuous Markov chain taking values in § = 
{1,2} with the generator T = (5^)2x2: 

-Qu = 9i2 = 1, -922 = 921 = q > 0. 

The unique stationary probability distribution of the Markov chain r(t) is 

7r=(7ri,7r 2 ) - (g/(l + «), 1/(1 + «)). 

Let /, g : R x S — > R be Lipschitz continuous in the first argument and satisfy linear 
growth conditions. Assume that there exist constants bi £ R, d% > 0, i = 1,2, such 
that for x £ R 



(3.2) 

and 

(3.3) 

For i = 1,2 let 

Si := 
m 4 := 

Assume further that 



2xf(x, i) + g 2 {x, i) < hx 2 , 
xg(x, i) > d? x 2 . 

[ 7l 2 (y) + 2 7l (y) - 21n(l + 7i (y))}X(dy) > 0, 

X) 

[2 7i (y)-21n(l+ 7i (y))]A(dy). 



(3.4) 



[(ln(l + 7 ,(y))) 2 + 7 2 (y)]X(dy) < 00, i = 1, 2. 



n 



Consider the following equation 

dX(t) = [AX(t) + f(X(t),r(t))} dt + g(X(t), r{t))dW{t) 



(3.5) +/ 1 (r(t),y)X(t-)N(dy,dt), t>0,x £ (O.tt); 

Jo 

X(0,ac) = zt (x), xe(0,7r); 0) = X(t, tt) = 0, t > 0. 
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In this example, set H := L 2 ([Q,ir]), A — ^( a ( a; )^) with domain V(A) satisfying 
the boundary conditions above. We let a(x) be a measurable function defined on 
[0, 7r] such that 



(3.6) 



< v < a(x), 



for some positive constant v. 

Let U(t,u,i) := Aj||u||^,u € H, i — 1,2, where Ai = 1 and A2 is a positive 

constant which will be determined later. Note from (13.21) . (|3.6p and Poincare's 
inequality that for u £ T>(A) 

CU(t, u, 1) = 2(Au + f(u, 1), u) H + \\g(u, 1)\\ 2 H 

+ / l 2 i(y)Hdy)\\u\\ 2 H + q n Xi\\uf H + gi 2 AaH||r 



< 



7 2 (y)X(dy) + \ 2 - I 



U{t,u,l) 



-2v + bi -\ 
:= axU(t,u, 1), 

and similarly 

CU(t,u,2)<[-2u + b 2 + J tf(y)\(dy) + q(J--i)]u(t,u,2) 
:= a 2 U(t,u,2). 
By the definition of U, it is easy to see that 



pi = 1 - A2 + In A 2 and p 2 = q(l — In A 2 

A 2 



From ([3~3")1 , it follows that 



BU(t,u,i) > ck\\u\\% = ^U 2 {t,u,i) := ^U 2 (t,u,i). 
A; 



Moreover, (v) follows from (|3.4I) and (iv), and (vi) holds due to the definition of U. 
Thus, by Theorem 13. II we arrive at 



^U V \H\\X{t)\\H)<-^y 



a.s. 



where 



d := q 
In particular, let 



di 



mi + 2^ — 61 



2X 2 



mi + 2v - bo 



y + mi + 2v - h < 0, 



and choose A9 > such that 



— 2 + m 2 + 2v - b 2 > 0. 
2A 2 

Then Eq. (|3.5p is almost surely exponentially stable whenever 

0<q<- (Jjp + m 2 + 2v - 6 2 ) I (y + Tin + 2^ - &! 
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4. Linear Switching-diffusion SPDEs with Jumps 



In this section to demonstrate that the results obtained in Theorem l3.1l are sharp, 
we will discuss a class of linear switching-diffusion SPDEs with jumps. 

For a bounded domain O C R" with C°° boundary dO, let H := l?{0) 
denote the family of all real-valued square integrable functions, equipped with 
the usual inner product (f,g)ii '■= J f(x)g{x)dx, f, g G H and norm := 

(Jo f 2 ( x )d x ) 2 j f £ H. Let A := be the classical Laplace operator de- 

fined on the Sobolev space V(A) := H%(0) n H 2 (0), where H m (0),m = 1,2, 
consist of functions of H whose derivatives D a u, in the sense of distributions, of 
order \a\ < m are in H. It is well known that there exists an orthonormal basis of 
H, {e„}„>i,n = 1,2,..., satisfying 

(4.1) e„ G V(A), -Ae n = X n e n . 

Thus, for any / G H, we can write 

oo 

/ = ^ f n e n , where /„ = (/, e„). 

Consider the following SPDE driven by a switching-diffusion process with jumps 

dX{t) = (AX(t) + a{r{t))X{t))dt + (3{r(t))X(t)dW(t) 
+ / °° j(r(t),y)X(t-)N(dt, dy), xeO,i>0, 
X(t,x) = 0,i£ dO, t>0, 
X(0, x) = u°(x), x G O and r(0) = r Q . 



(4.2 



Here a, f3 : S — > R, 7 : § x [0,oo) — > R, W is a real- valued Wiener process on the 
probability space {tt, T, P}. 

Applying Theorem 13. II we can deduce the following conclusion. 



Theorem 4.1. Assume that for i G S, y G (0, 00) 

poo 

(4.3) 7i (y) :=7(i,y) > -1, max / jf(y)X(dy) < 00, 
and 

/•OO 

(4.4) max / (ln(l + ^{y))f \{dy) < 00. 



Then Eq. ((4721) has the property 

limsup-ln(||X(t)|| H ) < -Ai+£>, + / Kl + 7<(l/)) -7<(l/)]A(dj/) 

In particular, the solution of Eq. (|4.2p is almost surely exponentially stable if 



Ar > J2 *j \&j - + J o [ln(l + 7i (y)) - UvMdy)) . 

Proof. Under (|4.3p . by Lemma |2~21 Eq . (|4.2I) admits a unique global mild solution. 
Let t/(u) := G H, and set for t > Q,u ^ H, i 6 § and yeZ 

.A(t, u, i) := Att + o^u, g(t,u,i) :— fiiU and $(£, u, i, y) := ji(y)u. 
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It is easy to compute cti , (3i and 5% in Theorem 13.11 as follows 

/>oo 

a. L = -2X l + 2a, + pj + / jf(y)X(dy), ft = 4ft 2 , 
Jo 

and 

/>OG 

Si= / [2(7i(tf)-Ml+7i(y)))+7?(»)]A(d»). 
Jo 

Moreover, noting that 

A(i, u, i, y) = (1 + li{y)) 2 and T(i, u, r , i, I) = 1, 
together with (|4.3I) and (|4.4I) . we can deduce that (iv) and (vi) hold. Then the 



desired assertion follows by Theorem 13. II □ 

We now discuss the sample path stability of the solution to Eq. (|4.2j) using its 
explicit mild solution, which will be given in Lemma l4.2l below. In the sequel, when 
u° is deterministic and u° ^ 0, we set no := inf {n : u„ ^ 0}, where it„ := (u , e n ) 
for n > 1. 

Lemma 4.2. Under (|4.3[) . the unique global mild solution of Eq. (I4.2[) has the 

explicit form 

(4.5) 

ft pOC 

X(tx) =v(tx)exvt -~ / /3*(r(s)Us+ / / IWl + 7fr(s), «)) - -yfrfs), «)lAWy)rfs 



:{t,x)=v(t,x)exp{-~J P 2 (r{s))ds + J J + j{r(s),y)) - j{r{s),y)]\{dy)ds 

+ f p(r(s))dW(s) + f f ln(l + j(r(s),y))N(ds,dy)}, 



where 

(4.6) v(t,x) := V]exp{-A n t + f a(r(s))ds\ u° n e n (x), t>0, xeO. 

n=l L ^ 



Proof. Under (|4.3j) . Eq. (|4.2|) has a unique global mild solution such that 
(4-7) 

X(t) = S(t)u°+ [ a(r(s))S(t-s)X(s))ds+ [ j3(r(s))S(t - s)X(s))dW(s) 



/ j(r(s),y)S(t- S )X(s-)N(ds,dy). 
o Jo 

Since 

oo 

S(t)u = e~ Xnt (u, e n )H£n for u <E H, 

n=l 

Eq. (I4.7P can be rewritten in the form 

ft 



X(t) = ^e- A "*( M °,e„)//e n + / a(r( S ))^ e - A "( t - s '(X( S ),e„) ff e„ds 

n=l n=l 

„t oc 

+ / ^(r( S ))^e- A "( t - s )(X( S ), e „) ff e„^( S ) 
"/O n=l 

OO 

7(r(s), y) £ e - A «(*" s ) (X ( s -), e n ) H e n N(ds, dy). 



n=l 

t />oo 



JO 
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This further yields 

e Xnt (X(t),e n )„ = (u°,e n ) H + f a(r(s))e x " s (X(s),e n ) H ds 

Jo 

+ ( P(r{s))e x - S (X(s),e n ) H dW{s) 
Jo 

rt poo 

7(r(s), y)e x " s (X(s-), e n ) H N(ds, dy). 



o Jo 



Then Ito's formula gives 



(X(t),e n ) H = (w°,e„) ff exp{ - A„i + ^ [a(r(s)) - ^ 2 (r(s)) 

t (• do 

[ln(l + 7(r(s), y)) - j(r(s),y)]X(dy)ds 



o Jo 

t rt poo 



+ P(r(s))dW(s)+ / ln(l + 7(r(s),y))N(ds,dy)\. 

Jn Jo Jo > 

He follows by observing tl 

oo 

x(t) = J2(x(t), 



>0 JO 

The desired assertion therefore follows by observing that 



□ 



Theorem 4.3. Under the conditions of Theorem 14.11 Eq. (14.21) has the following 
properties: 



(i) If u° is deterministic and u° ^ 0, then 

i m / | 

hm - ln(\\X(t)\\ H ) = -X no «i - 'J. 



2 

+ / [ln(l + li{y)) - Ji(y)]X(dy) ) a.s. 



In particular, the solution of Eq. (|4.2I) with initial condition u° will converge 
exponentially to zero with probability one if and only if 

m / 1 - f°° \ 
A„ > J2 *i If* ~ 2^ + J [ln(1 + ^ V)) ~ iMMW) ■ 

(ii) For any initial condition u , 

1 m / 1 - 

limsup-ln(||X(t)||K) < -Ai + ^} 

3 = 1 x 

+ / [ln(l + 7<(»)) - li{v)]Kdy) ) a-s- 



In particular, the solution of Eq. (|4.2[) is almost surely exponentially stable 
if 



Ai >E^ («J- + X [ln(1 



■ 7i(y)) - 7»(y)]A(dj/) 
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Proof. Note that 



±M\\X(t)\\ H ) = jH\\v(t)\\ H ) - I f 'p 2 (r(s))ds 



1 



jMx(t) + jM2(t), 



[ln(l + 7(r(s), y)) - j(r(s), y)]X(dy)ds 



where 



ft />oo 

Mi(t) := I P(r(s))dW(s) and M 2 (t) := I / ln(l + 7 (r(s), j/))JV(ds, dj/). 



Since 



o Jo 



(Mi) t < i max pf and (M 2 ) t < t max / (ln(l + 7*(v))) A(dy), 

KKm l<i<m 



together with the strong law of large numbers for local martingales, e.g., Lipster 
[IS] , we have 

lim -Mi(t) = lim -M 2 (t) = a.s. 

t— »oo £ t— foo £ 

Furthermore, by the ergodic property of Markov chains, e.g., |20[ Theorem 3.8.1, 
pi 26] 



i / ^( r ( s ))ds + ] f r[Ml+j(r(s),y))-*f(r(s),y)}\(dy)d 



2t L t 



o Jo 



5>(- 5 *+j£ Wi 



+ 7i(y)) -7i(y)]A(dy) ) a.s 



whenever t — >■ oo. Moreover, it is not difhcult to show that 

1 . ... ,. , J = -A„ + 537^1 ^i^' if ""° is deterministic , 
lim - ln(||i;(t)||#)< ^ 3 _ 

t->oo f I < Ai + 2_/j=i 7T j a jj f° r an y initial condition u . 

The proof is therefore complete. □ 

Example 4.4. Let r(t) be a right-continuous Markov chain taking values in § = 
{1,2} with the generator T = (3,^)2x2 : 

-?n = 912 = v > 0, -g 2 2 = 921 = q > 0. 

Then the unique stationary probability distribution of the Markov chain r(t) is 

7T = (tti, tt 2 ) = + g), + q)). 

For j = 1 , 2 set 

fii := / [ln(l+7i(y)) - ii(y)]\(dy). 



Let ai, «2 £K such that 

(4.8) ai +m > 1, a 2 + M2 < 1, 
and choose q obeying 

(4.9) < q < i/(l - a 2 - /J, 2 )/(ai + m - 1). 
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Consider the switching-diffusion equation with jumps: 

dX{t) = (AX(t) + a(r(t))X(t))dt 

+ / °° j(r(t),y)X(t-)N(dt, dy), t > 0, x G (0, tt), 

x(t,a) = x{t,w) = o, t>o, 

X(0,x) = u°(x) = x G (0,tt). 

The previous stochastic system can be regarded as the combination of two equations 

/>oo 

(4.11) dX(t) = (AX(t) + a 1 X(t))dt + / jx(y)X(t-)N(dt,dy), t > 0, x G (0, tt), 

Jo 

and 

/>oo 

(4.12) dX(t) = (AX(t) + a 2 X(t))dt + / */ 2 {y)X{t-)N(dt, dy), t>0,x£ (0, tt) 

Jo 

with the same Dirichlet boundary condition and initial condition as Eq. (|4.10l) . 
switching from one to the other according to the law of the Markov chain. Recall- 
ing that e n (x) — *\/2/7rsin nx, n = 1, 2, 3, ... , are the eigenfunctions of —A, with 
positive and increasing eigenvalues X n = n 2 , we hence have Ai = 1. By Theorem 
14. li we have 

limsup - ln(||Jf (t)||#) < — 1 + tt\ol\ + 7r 2 a 2 + 7Ti/ii + 7r 2 /i2 a.s. 

t— >OG t 

This, together with (|4.8[) and (I4.9p . yields that Eq. (14. 10)) is almost surely exponen- 
tially stable. On the other hand, note that the initial condition u°(x) = -\/2/7rsina; 
is deterministic and v,® = 1, which implies no = 1. By Theorem 14.31 the solution 
to Eq. (|4.11[) has the property 

lim - ln(||X(t)||fl-) = —1 + Q.X + fix, a.s., 

t— yco t 

and the solution to Eq. (|4.12|) has the property 

lim -ln(||X(t)||jj) = -1 + a 2 + A*2, a -S- 

t— too t 

Then, by (|4.8p the solution of stochastic system (|4.11|) explodes exponentially, and 
the solution of stochastic system (|4.12p converges exponentially to zero. However, 
we observe that due to the Markovian switching the overall system (I4.10|) is almost 
surely exponentially stable. 
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